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We have used a combination of numerical modeling and experiments to study carbon etching in 
the presence of a hydrogen plasma. We model the evolution of a low density EUV-induced plasma 
during and after the EUV pulse to obtain the energy resolved ion fluxes from the plasma to the 
surface. By relating the computed ion fluxes to the experimentally observed etching rate at various 
pressures and ion energies, we show that at low pressure and energy, carbon etching is due to 
chemical sputtering, while at high pressure and energy a reactive ion etching process is likely to 
dominate. 


I. INTRODUCTION 

Many vacuum mirror optical tools suffer from the 
build-up of carbon contamination due to cracking of 
hydrocarbons under powerful vacuum ultraviolet radia¬ 
tion [1]. The large absorption of EUV radiation by car¬ 
bon becomes significant in the case of multi-element op¬ 
tical systems, where throughput can be greatly reduced 
by even a very thin layer of carbon contamination on the 
top of each mirror. 

The problem of EUV induced carbon contamination 
has been addressed in a series of publications [2-5] . The 
reduction of carbon films in a hydrogen atmosphere or hy¬ 
drogen plasma has also been extensively studied [6-10]. 
Despite numerous studies, however, it is still difficult to 
predict the carbon removal rate, because there are many 
contributing factors. Several aspects that significantly 
affect the carbon removal rate include: many different 
allotropes and compounds of the carbon (e.g. soft black 
or hard graphite), many different contributing reaction 
paths (e.g. physical sputtering, chemical sputtering, re¬ 
active ion etching etc.), and, last but not least, small 
admixtures to the background gas, which, while residing 
on the carbon surface, can produce reactive species once 
irradiated by EUV. Despite this complexity, experiments 
have shown that carbon etching can be achieved under 
certain EUV-induced plasma conditions. Nevertheless, it 
has proven difficult to fully understand the etch process, 
because the characteristics of the EUV-induced plasma 
are poorly known, and the plasma-surface interaction has 
many contributing factors [10]. 

In this paper, we use a model of the EUV-induced 
plasma to numerically analyse the fluxes from the plasma 
to the sample surface. Our model is a self-consistent 2D 
Particle-in-Cell model of the weakly ionized low pressure 
hydrogen plasma that is formed during the EUV pulse 
due to ionization by EUV photons and secondary elec¬ 
trons from the surface. As described in Section III, our 
model provides an accurate estimate the ion flux compo¬ 


sition and energy distribution. However, a considerable 
number of parameters required for accurate simulations 
are not well known, therefore, the modeling results were 
combined with experimental observations. With the aid 
of simulations, we show that the shape of the energy dis¬ 
tribution function of the ion fluxes in the considered ex¬ 
perimental setup are mainly defined by the setup geom¬ 
etry, background pressure, and externally applied bias 
voltage. The ion dose, on the other hand is sensitive to 
the variations of many other parameters (e.g. EUV dose, 
secondary electron yield, etc.). 

By combining the computed energy distribution func¬ 
tion of the ion flux with the experimentally measured 
ion dose, insight into the mechanism for carbon removal 
was gained. By analysing the differences in yield be¬ 
tween EUV-induced plasma and surface wave discharge 
plasma experiments in combination with numerical sim¬ 
ulations, we show that chemical sputtering dominates for 
low pressures and energies. It was found that the carbon 
removal yield for both the surface wave discharge and 
EUV-induced plasmas was similar in the overlapping en¬ 
ergy range. Hence, the effect of the EUV radiation on 
carbon removal is found to be significantly smaller than 
was estimated previously [10]. 

II. EXPERIMENTAL SETUP 

The ISAN EUV experimental setup is based on a 
tin EUV radiation source, which is a Z-pinch discharge 
plasma with 1500 Hz repetition rate, which has been de¬ 
scribed in detail elsewhere [11]. This source is a good tool 
for exploring EUV-induced surface processes over a large 
number of pulses (> 1 MShot). In brief, EUV radiation is 
introduced into a so-called ’’clean” chamber (see Fig. 1), 
separated from the source and collector optics by a Si:Zr 
spectral purity filter (SPF). The clean chamber is differ¬ 
entially pumped to pressure of 3.5 x 10 -8 torr. Under 
vacuum conditions, the background hydrocarbon carbon 
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Figure 1: Configuration of the experimental chamber. The 
metal cylinder inner diameter is 28 mm, the sample holder 
diameter is 18 mm, the inner radius of dielectric mica di¬ 
aphragm is 8 mm, the distance from sample to grid is 2.4 cm, 
distance between the grid and SPF is 1.5 cm. 


growth rate was measured to be 0.4±0.2 nm/10 Mshot. 

The diameter of the EUV beam at the sample was 
5 mm. In addition to the direct beam, some scattered 
EUV radiation was also incident on the sample. The 
EUV pulse duration is about 100 ns (FWHM), with a 
tail, as described in [12]. 

The incident EUV power was measured using a sen¬ 
sitive thermo-couple attached to a thin copper disk. It 
was found that the radiation intensity was 0.13 W/cm 2 
after the SPF and approximately 0.75 W/cm 2 without 
the SPF. The ratio of EUV intensities with and without 
the SPF corresponds to the calculated transmission of 
the 100 nm Si:Zr SPF filter over an energy range of 60 
to 100 eV (see Fig 2). 

The EUV intensity on the sample decreases with time 
because of carbon growth on the SPF filter and focus¬ 
ing optics. In later experiments, the EUV intensity was 
measured to be 0.1 W/cm 2 after the SPF. 

During the experiments, the hydrogen pressure was set 
in the range of 2.8 Pa - 86.5 Pa. To ensure that each 
radiation pulse excited a plasma in an atmosphere dom¬ 
inated by hydrogen, hydrogen flowed through a liquid 
nitrogen trap (to remove water) and into the chamber at 
100 liters xtorr/minute . 

To control the energy of the ion flux, the sample holder 
assembly was biased in the range of -200 - 0 V, while all 
other metallic electrodes were grounded. The samples 
consisted of carbon, deposited by magnetron sputtering 
on a silicon wafer to a thickness of ~30 nm. Each sample 
was exposed to the 10' EUV pulses at a different com¬ 
bination of bias and hydrogen pressure. After exposure, 
the amount of carbon removed was measured by X-Ray 
fluorescence (XRF) (EDS) analysis and spectroscopic el- 
lipsometry. 


III. MODEL 

A two dimensional particle in cell (PIC) model with rz 
geometry was used to model the experiment. Our model 
follows the general PIC scheme, described elsewhere [13]. 

Ionization is initiated by an EUV pulse, which directly 
ionizes the background gas, and produces electrons by 
photo emission from the SPF and the sample. This pro¬ 
cess leads to the formation of an EUV induced plasma in 
the chamber. The ionization and photo-emission process 
and their inclusion in the model is described below. 

Although the plasma is continually re-ignited by the 
pulsed EUV radiation, we consider each EUV pulse to 
induce a plasma in a cold neutral gas that is at equi¬ 
librium, because the characteristic time of plasma decay 
in the given geometry is about ~ 20/zs, which is much 
shorter than the 660 /zs between pulses. The restora¬ 
tion of thermal equilibrium in the background gas is also 
much shorter then the time between pulses. 

A. Chamber configuration 

The configuration of vacuum chamber is presented in 
Fig.l. Because the chamber is axially symmetric, it is 
possible to include the structure of the internal chamber 
in the simulations. To accurately model the plasma dy¬ 
namics when the sample is biased, the space between the 
SPF and the metal grid is included in our simulations. 
The grid is included in the model as a number of metallic 
rings. When a bias is applied to the sample, ions from 
the plasma, which is formed between SPF and metal grid, 
are pulled towards the sample. 

B. Deielectric model 

In the experiment, a dielectric mica diaphragm was 
used to prevent EUV radiation being incident on the sam¬ 
ple holder structure, thereby decreasing electron emission 
from the sample holder structure. Therefore, allowing 
plasma parameters to be estimated from the discharge 
characteristics of the sample. In the model, this feature 
is included with a simple dielectric model: a dielectric 
is not conductive, it can accumulate charge and the sec¬ 
ondary electron yield under EUV irradiation is order of 
magnitude lower than for the sample. 

The electron reflection and secondary electron emis¬ 
sion under electron impact (SEE) are included in the 
model. These two processes are combined into SEE with 
the probability ( P(E )) defined as follows. 

P(E) =p 0 + ^(1.0 -po) (1) 

Here E - is the energy of incoming electron. The param¬ 
eters of the used mica are unknown, therefore, we choses 
Ei = 45 eV and p 0 = 0.7, these parameters corresponds 
to linear fit for SEE yield for quarts in [14]. 
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This approximation blends between significant slow 
electron backscattering for low energy electrons and sec¬ 
ondary electron emission for higher impact energies. We 
do not use higher order approximations, since both ef¬ 
fects are known to depend strongly on the surface condi¬ 
tions [14, 15]. 

This model for the dielectric is over simplified, but suf¬ 
ficient to allow us to to take the out of focus EUV radi¬ 
ation, which is incident on the diaphragm, into account 
in a consistent manner. 


Here P(E) is the probability of emitting an electron with 
energy, E, from a surface with a work function, IV. For 
carbon TV = 5 eV [19]. For electrons emitted from the 
surface, we assume an angular dependence given by a 
cosine emission law [18]. For purpose of estimation, we 
used the same approximation for the spectrum of photo¬ 
electrons emitted from the dielectric. 


E. EUV spectrum and photoionization 


C. Length scales and grid resolution 


The time evolution of one EUV pulse has two distinct 
stages: negative space-charge dominated during the be¬ 
ginning of the EUV pulse, and decay of the positively 
charged plasma after the pulse [16]. 

In order to model the negative space-charge dominated 
part of the plasma evolution, the potential well near the 
surface, which is formed by the electrons that have es¬ 
caped to the volume, must be resolved. The length scale 
of the space-charge potential well depends on the energy 
distribution of photoelectrons from the surface and on 
the current which passes through the system. 

For the purpose of estimation, it is possible to simplify 
the formulas from [17] and obtain 


z m ~ 0.1 cm x 


(T[eV]) 3 / 4 

(/[mA/crn 2 ]) 1 / 2 


( 2 ) 


Here, T is the initial temperature of the emitted pho¬ 
toelectrons, / is the current density near the cathode, 
and z m is the distance from the cathode to the bottom 
of the space charge potential well. For T ~1 eV and 
I ~20 mA/crn 2 , one obtains z m ^0.02 cm. 

For the plasma dominated part of the plasma evolu¬ 
tion we need to resolve the Debye length in the volume, 
and also have a finer grid near the surface to resolve the 
plasma sheath. 

From the EUV intensity (~0.02 mJ per pulse) and 
background hydrogen pressure (2.8 Pa, 11.2 Pa), the 
plasma density can be estimated to be about 10 9 cm -3 . 
For T e ~0.5 eV and N e ~ 10 9 1/cm 3 one obtains 
Rb ~0.015 cm, which is comparable to the space charge 
length scale estimated above. To achieve this, the grid is 
refined near the sample surface and SPF. The minimum 
grid cell size is 0.05 x Rb, and is gradually increased in 
steps of 5% until the bulk cell size of 0.5 x R D is reached. 
Several tests were performed to ensure that the chosen 
grid resolution does not affect the plasma dynamics. 


To calculate the distribution of ion species, the spec¬ 
trum of the EUV radiation must be included in the 
model, along with the energy-dependent cross sections 
for each photoionization process. Direct photoionisation 
of hydrogen by EUV photons is included as two processes: 

H 2 + hu —y H[[” + e 
H 2 +/u/ -> H++H + e 

The photoionisation cross-sections were taken from [20]. 

In the experiment, the EUV spectrum was measured 
(green curve Fig. 2) and found to be in agreement with 
data from [12]. The EUV spectrum, however, was mea¬ 
sured before reflection from the Zr collector mirrors, and 
transmission through the SPF. 

The calculated transparency curve that was provided 
with the Si:Zr SPF was used to calculate the EUV spec¬ 
trum after transmission through the SPF (See Fig. 2 blue 
curve). 

The SPF transmission curve has a transmission of 
about 1% in the range of 20 30 eV, which is important 

because the photoabsorption cross-section is very large 
in this energy range (see Fig. 2 red curve). 

The contribution due to radiation in the 20 - 30 eV 
range is difficult to quantify for two reasons: the inten¬ 
sity of the radiation varies significantly due to carbon 
growth on the SPF, and the accuracy of the SPF trans¬ 
mission curve in this energy range is unknown. Never¬ 
theless, we include it because, even a small transmission 
will result in additional photoinization, which leads to a 
reduced space-charge potential barrier at low pressures 
(e.g. 2.8 Pa, 11.2 Pa). 

Although the photoionization due to the 20 - 30 eV ra¬ 
diation range is significant compared to the 13 nrn band, 
it has significance only for cases where no bias is applied 
to the sample. For cases with bias, the presence of this 
radiation leads to a decrease of the space-charge poten¬ 
tial. 


D. photo-electron emission 


F. Cross-sections set 


The energy spectrum of photoelectrons emitted from 
the SPF and sample surface are calculated from [18]: 


P(E) 


E 


( 3 ) 


The cross-sections set used in the model consists of 
electron collisions with hydrogen and ion collisions with 
hydrogen. The collisions between plasma species and 
three body processes are neglected due to their low prob¬ 
ability under the conditions considered here. To accu- 


(E + TV) 4 
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Figure 2: EUV spectrum after SPF as used in the simulations. 
Note, that the small SPF transmission, estimated at ~1% of 
EUV energy per pulse, significantly contributes to the number 
of photoionisation events in the volume. 



Bias (V) 


simulation 

Ar—A 0 Pa 0—<) 22.4 Pa 

2.8 Pa ►—►45 Pa 
*-★ 11.2 Pa 


experiment 
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O O 2.8 Pa >- > 45 Pa 
*- * 11.2 Pa 


Figure 3: Collected charge as function of bias, 0.017mJ/pulse 
(0.13 W/cm 2 at 1500 Hz to the spot of 5 mm diameter. 


rately model electron collision related processes in hy¬ 
drogen discharges with a Monte Carlo (MC) model, one 
needs to take into account the differential cross-sections 
for ionisation and excitations processes. As described in 
detail in [21], the particular choice of the angular depen¬ 
dence of cross-sections significantly influences the simu¬ 
lation results. We adopt a set of cross sections found in 
Ref. [21]. 

We use an approach similar to that described in [21]. 
The reaction probability is sampled using the integrated 
reaction cross-section. The differential cross-section data 
is used to determine the collision kinematics and energy 
redistribution between products. Electron elastic scat¬ 
tering and hydrogen electronic excitations, and angular 
scattering data is taken from [22]. For electron impact 
ionization of hydrogen, we use an experimentally deter¬ 
mined doubly differential cross-section [23, 24]. The set 
of cross-sections for collisions between ions and hydrogen 
is based on [25] because this set provides good agreement 
with swarm data for ions in hydrogen. We neglect the 
formation of H _ , because the cross-section of dissocia¬ 
tive electron attachment is very low, and the density of 
vibrationally exited hydrogen molecules too low to make 
a significant contribution to the production of H~. 

We make use of the procedure described in [26] to per¬ 
form Monte Carlo collisions with the background gas. We 
tested the consistency of our implementation by model¬ 
ing swarm experiments and found good agreement with 
experimental values [27] for the first Townsend electron 
ionization coefficient, the electron mobility, and for H + 
and H[j~ mobility in hydrogen [28]. 


IV. ANALYSIS OF THE CHARGE — BIAS 
CHARACTERISTIC 

Let us begin with the analysis of the charge — bias 
voltage characteristics (Q - V), which were measured 
during the experiments. The ion dose on the sample per 
pulse for the range of pressures and bias voltages consid¬ 
ered here was estimated from the experimental results as 
a difference between the collected charge for given pres¬ 
sure and bias voltage and the collected charge for vacuum 
conditions [10]. The measured curves are presented in 
Fig. 3. It is worth noting that these Q — V characteristics 
have some peculiarities. 

Firstly, all the characteristics reach saturation, and the 
turning point for the almost all curves is approximately - 
75 V. If there is a significant contribution from ionization 
by the accelerated secondary electrons from the sample in 
the drift regime, one would expect a rapid increase of the 
collected charge for biases in the range of -200 .. -50 V. 

Secondly, the increase of the collected charge for -200 V 
bias does not depend linearly on pressure in the range of 
2.8 Pa to 22.4 Pa, i.e. the collected charge increases less 
than two-fold while the pressure increases eight times. 
But, one would expect at least linear growth of the 
collected charge with pressure if the volume ionization 
provides a significant contribution to the total collected 
charge. Thus, the overall contribution of photoionization 
in the volume is small compared with the other factors. 

Thirdly, for -200 V bias and 22.4 Pa background pres¬ 
sure, due to the photoionization process, all the external 
bias would be applied over a very small layer near the 
sample, thus, all the secondary electrons will gain 200 eV 
energy. The H 2 ionization cross-section for this energy is 
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approximately 6 • 10~ 1 ' cm 2 , thus, in the space between 
the sample and the grid (l ~ 2.5 cm), these electrons 
will, on average, have 0.8 ionization inducing collisions. 
If we assume that, under vacuum conditions and -200 V 
bias, the SE current from the sample saturates, and we 
obtain that, just due to the direct electron induced ion¬ 
ization, without cascade process and photoionization, the 
collected charge should be greater than 1.8 nC. in the ex¬ 
periment, however, the collected charge is approximately 
1.9 nC. 


A. Average secondary electron yield 

The above analysis suggests that the combination of 
EUV power per pulse, and effective secondary electron 
yield (SEY) produces approximately 1 nC of electrons 
from the sample. 

It is instructive to estimate the total SE charge from 
the experimental parameters. From the EUV intensity, 
repetition rate and spot size, the average dose per pulse 
was approximately ^0.017 mJ. The secondary electron 
yield for carbon under EUV radiation is estimated to be 
approximately 'ysE ~ 0.01 for the photon energy range of 
60 eV - 100 eV [2, 29], which leads to a saturation charge 
of 2 nC, which would lead to a significant disagreement 
between the simulated and measured charge bias charac¬ 
teristics. It is also worth noting that effective SEY from 
mica was small, since, in the experiment, the charge col¬ 
lected from a sample made from mica was measured to 
be an order of magnitude smaller than for the carbon 
sample. 

However, the measured charge bias characteristics are 
the only experimental data that provides a reference 
point for the simulations of the experiment dynamics. 
Since the spatial distribution of EUV intensity is sub¬ 
ject to an unknown systematic measurement error, and 
the SEY is known to vary widely, depending on the sur¬ 
face conditions, we chose to keep the product of incident 
EUV and SEY a constant, chosen to provide 1 nC total 
SE charge from the sample. 


B. Role of dielectric ring 

During the simulations we found that the mica accu¬ 
mulated charge over many EUV pulses, significantly ef¬ 
fecting the local field distribution, and, hence changing 
the flux incident on the sample. If charging is neglected, 
the Q - V response under vacuum and 2.8 Pa conditions 
cannot be reproduced for any reasonable parameter val¬ 
ues. The discrepancy is caused by the potential barrier 
near the sample surface, which is created by the negative 
space charge generated during the EUV pulse. Although 
the potential is low, it produces a significant effect be¬ 
cause the SE energy spectrum (3) is strongly peaked at 
rather low energies: approximately W /3 = 1.3 eV. 



R (cm) 


Figure 4: Comparison of the potential distibution hear the 
sample holder for uncharged mica (left) and pre-charged mica 
for -100 V bias on the sample holder. The charge on the mica 
significantly increases the local field strength near the smaple 
surface. 


The charge accumulation on the mica significantly in¬ 
creases the local field strength near the sample surface 
(see Fig. 4). This removes the space charge potential 
barrier for biases of -100V bias and higher under the con¬ 
ditions that we consider. 

The charging of the mica should saturate, leading to an 
unchanging charge density distribution on the mica. To 
estimate the charge density, we simulated one hundred 
pulses under vacuum conditions for all experimentally 
applied bias values. In these calculations, the effective 
SEY from the mica is assumed to be about 0.001, which 
corresponds to the experimentally measured value. 

This approach corresponds to the experimental proce¬ 
dure, since the Q - V characteristics were measured via 
averaging a large number of pulses for every combination 
of bias and pressure. No special means were used in the 
experiment to remove the accumulated charge from the 
mica between pulses. 

For negative biases, simulations show that the mica po¬ 
tential rapidly reaches approximately 0 V potential. For 
these conditions, further charging is very slow, because 
(a) the electric field strength near the mica is very small, 
therefore, even a small amount of emitted charge creates 
a space charge potential barrier, and, (b), because the 
mica potential is close to zero, there is a small flux of 
electrons from the sample to the mica. 

The continued slow variation in charge distribution is 
impractical to simulate, however.The final charge distri- 
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bution requires a very long time to calculate and is sensi¬ 
tive to the combination of the SEY from both the sample 
and the mica, and to the corresponding SE energy distri¬ 
bution functions. Nevertheless, the additional errors due 
to these limitations are expected to be small. In practice 
the maximum potential of the mica under vacuum condi¬ 
tions is bounded, since, as the mica’s potential increases 
to ~10 V, the electron current from the sample to mica 
becomes significant. 

It is worth noting that the EUV plasma itself can 
contribute to charging the mica, since, for stationary 
discharges, the dielectric charges to the plasma poten¬ 
tial. Despite the the fact that the potential of the EUV- 
induced plasma can be high (e.g. 20 - 30 V or more after 
EUV pulse), it rapidly decreases to several volts due to 
electron cooling due to collisions with the background 
gas. Most probably, the mica potential would converge 
to the time averaged potential of the plasma, e.g. several 
volts. Therefore, in the following results, the mica was 
pre-charged to the value that was obtained from simula¬ 
tions under vacuum conditions for the appropriate bias. 

The optimum parameter combination that reproduces 
the experimentally measured Q - V characteristics were 
chosen as follows: the mica was pre-charged as described 
above, the EUV intensity was kept constant at the mea¬ 
sured value, but the amount of the scattered EUV radi¬ 
ation incident on the sample holder was chosen so that 
40% of EUV radiation was incident on the dielectric mica, 
and the effective SEY from the sample was kept at to 0.01 
for the photon energy range of 60 eV - 100 eV, while 
the SE contribution from the VUV part of the spectrum 
was neglected. The simulated Q - V characteristic, after 
parameter optimization, and comparison with the mea¬ 
surements are presented in Fig. 3. 


V. ION FLUXES TO THE SAMPLE SURFACE 

The removal of carbon from the surface material should 
be directly dependent on the energy distribution function 
(EDF) of the ion flux incident on the surface. The con¬ 
ditions of the experimental study of carbon cleaning in 
the ISAN EUV experiment are summarized in Table I. 
We computed the plasma conditions, ion fluxes, and ion 
EDF for these conditions. The simulated time for all 
cases was 10 /zs. For 2.8 Pa - 11.2 Pa and biases -200 V 
- -100 V, the plasma had completely decayed and all ions 
were collected during the simulation time. But, for the 
60 Pa and 86.5 Pa cases some plasma was still left in the 
simulation domain. 

The EDFs, integrated over the simulation window, and 
averaged over the sample surface, are presented in Fig. 5. 
The composition of the computed ion flux depends signif¬ 
icantly on the pressure. For the 3 Pa case, the main ion is 
H 2 , because the characteristic time for to H(j~ conver¬ 
sion in 3 Pa H 2 is approximately 0.5 fi s. With increasing 
pressure, H(|" becomes the main ion, as expected. There 
is also a non-negligible contribution from fast H 2 which 


is produced due to a resonant charge exchange reaction. 

The maximum ion energy is larger than the applied 
bias voltage due to the build-up of the plasma potential. 
But the number of such fast ions is small, due to the fast 
plasma decay. For the same reason, the most energetic 
ion is H + because its small mass allows it to accelerate 
during the decay of the plasma potential. 

It was observed, that contrary to the ion dose, the 
shape of the ion flux EDF was mainly defined by bias, 
pressure and the accumulated charge on the dielectric. 
However, as discussed in previous section, we expect that 
the equilibrium mica potential does not significantly dif¬ 
fer from 0 V. For small variations of the mica potential (if 
all other parameters are the same), the variations of the 
EDF shape are small for all simulations where a sample 
bias was applied. Moreover, the EDF was barely sensitive 
to reasonable variations of other parameters. 

This insensitivity is a consequence of the experimental 
procedure. The applied bias allows ions to be collected 
from the entire chamber volume. Hence, the characteris¬ 
tics of the plasma are important only for the short time 
when the sheath between plasma and sample is small. 
As ions are collected at the sample, the sheath size be¬ 
comes larger and larger. Therefore, as time progresses, 
the instantantaneous ion flux EDF starts to depend only 
on the bias potential, mica potential, and the distance 
that the ions travel between the plasma and the sample. 
This is because these parameters determine the maxi¬ 
mum possible energy of the ions, while the background 
pressure determines the number of collisions (e.g energy 
loss), and, hence, the average energy of the ion at the 
sample. Therefore, for the considered conditions, the 
shape of the ion flux EDF is mostly determined by the 
bias, mica potential, pressure, and the chamber geome¬ 
try. 

It is notable, that the mica charge, due to exposure to 
EUV radiation and plasma, leads to the ion flux being 
focused on the sample (see Fig. 7). The ion flux to the 
mica itself becomes negligible compared to the simulation 
with zero charge on the mica. Therefore, the measured 
ion doses should represent accurately the ion doses on 
the sample. 

Hence, for the analysis of the experiments, it is rea¬ 
sonable to take the computed EDF shape, normed to the 
experimentally measured ion dose (if it is available). This 
approach allows us to compare the experimental results 
at different pressures. 


VI. DISCUSSION 

The maximum etched depths and total etched volumes 
per exposure of 10 7 pulses for different pressures and bi¬ 
ases are presented in Table. I. The etched volumes are 
the integrals over etch profiles, which were determined 
by XRF measurement (see Fig. 8). In the cases of 60 Pa 
and 86.5 Pa pressure, the etch profile was only measured 
at the center of the EUV spot. We assumed, for these 
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Figure 5: Energy distribution function of ion and fast neu¬ 
trals flux to the surface for 3 Pa (top) and 11.2 Pa (bottom) 
for -100V bias. The fluxes were integrated over time and over 
sample surface and sum of fluxes normalized to unity. Note 
that Hg" spike is on the same position over energy in spite of 
the increase of pressure. 

two cases, a uniform etching profile to estimate the car¬ 
bon removal yield. Hence, for these cases, the amount of 
removed carbon is most probably overestimated. 

As was discussed in the previous section, we used the 
experimentally measured ion doses for yield estimation. 
However, for some experiments, the ion dose was not 
measured. In these cases, we used the ion doses obtained 
from simulation results. This approach is valid at high 
pressures because the uncertainty in the simulations is 
on the order of 30% (see Fig. 3). 

For all samples, the average carbon removal yield was 
larger than 0.2 carbon atoms per ion. Thus, physical 


Figure 6: Comparison of the shape of the ion flux EDF ob¬ 
tained with pre-charged dielectric mica for different pressures. 
The ion fluxes were integrated over time and over sample sur¬ 
face, summed over ion types and normalized to unity. The 
carbon removal yield curve corresponds to surface wave dis¬ 
charge plasma (SWD) from [10]. The dotted part of the yield 
curve corresponds to linear extrapolation. 
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Figure 7: Comparison of the ion flux spatial distribution ob¬ 
tained with pre-charged mica and with zero charge on the 
mica. The ion fluxes were integrated over time and energy. 
Both distribution are normed to one if integrated over sample 
and dielectric mica surface areas. The conditions are 3 Pa 
background pressure and -100 V bias on the sample. The 
charge on the mica focuses the ion flux on the sample. Other 
combinations of pressure and bias show similar behavior. 


EUV intensity, mJ/(cm 2 x pulse) 













































































































Table I: Experimental conditions and carbon removal rate 


pressure 

Pa 

bias 

V 

EUV power 
W/cm 2 

removed carbon 

max etched depth etched volume 

nm/10 7 pulse 10 -6 cm 3 /10'pulse 

ion dose 6) 

nC/pulse 

average yield 

present study SWD recomputed 
C atom/ion C atom/ion 

3.0 

-50 

0.11 

1.7 ± 0.5 

~ 0.03 

~ 0.1 

0.5 

0.4 

3.0 

-100 

0.11 

5.7 ±0.5 

~ 0.1 

~ 0.25 

0.6 

0.5 

3.0 

-200 

0.13 

14 ±0.5 

~ 0.3 

~ 0.25 

1.9 

— 

11.2 

-200 

0.13 

30 a) ±0.5 

~ 0.5 

~ 0.5 

>1.6 

— 

60 

-100 

0.13 

23 ±0.5 

~ 1.2 c) 

> 2.2 d) 

<0.9 

0.5 

60 

-200 

0.1 

27 ± 0.5 e) 

~ 1.4 c) 

> 2.3 d) 

~1.0 e) 

— 

86.5 

-100 

0.2 

19 ±0.5 

~ l c) 

> 3.0 d) 

<0.5 

0.4 


a) The substrate was exposed, thus only lower boundary for cleaning rate could be estimated. 

b) Estimated from measured Q - V characteristics (see Fig. 3) as the difference between the collected charge for a given 
pressure and bias voltage, and the collected charge for vacuum conditions. 

c) The carbon removal was measured only in the center of the exposed area. However, an SEM image of the EUV exposed 
area showed no features. Therefore, for the purposes of estimating the carbon removal rate, we assume that the carbon was 
etched uniformly over the exposed area of the sample. 

d) Experimental Q - V values were missing, therefore the simulated ion dose were used for estimation. Only ions which hit 
the sample during the simulated period are included. No corrections were made for the plasma which was present in the 
chamber after the end of the simulation. 

e) This value is only an order of magnitude estimation, because all carbon was removed from the sample. 


sputtering is not the main process, since for energies be¬ 
low 200 eV, the expected sputtering yield is lower than 
5 • 10 -2 carbon atoms per hydrogen ion [30]. Therefore, 
carbon is removed via a chemical sputtering process, or 
a reactive ion etching process, because the effective yield 
for these processes is known to be large under certain 
conditions [6]. 

However, the results for -200 V bias, especially at low 
pressure (2.8 - 11.2 Pa), cannot be described by a chem¬ 
ical sputtering process, since the average yield is larger 
than one carbon atom per ion. Hence, reactive ion etch¬ 
ing may be responsible for these very high etch rates. 
However, for reactive ion etching there should be weakly 
bound radicals on the top of the carbon, which are des¬ 
orbed from surface due to ion impact, which, in turn, 
lead to a very high yield per ion. Although, we have no 
data about these radicals, it is likely to be some form of 
methane radical (e.g., CH^). The alternative: oxidation 
(e.g., CO and CO 2 ) [2] is unlikely because water, which 
is the main source of oxygen, is removed from the hy¬ 
drogen flow, and, during experiments, the chamber walls 
are kept at liquid nitrogen temperatures, trapping the 
majority of the residual water. Hence, the contribution 
of water to the carbon removal was limited by the initial 
coverage which is typically sub mono layer for amorphous 
carbon for the considered temperature range [31]. 

In the case of experiments with -50 V and -100 V bi¬ 
ases, the results can be compared with carbon etching 
in a surface wave discharge plasma [10], where similar 
magnetron deposited carbon on silicon substrate sam¬ 
ples were used. The carbon removal yield for experiments 
with biases of -50 V and -100 V was estimated by con¬ 


volving the simulated EDF with the energy dependent 
yield, as measured in the SWD experiments (see Fig. 6 
and column ”SWD recomputed” in table I). In spite of 
H]j~ being the main ion in SWD plasma, we used the same 
yield for estimation, since we do not expect dramatic dif¬ 
ference between Hij" and Hj for carbon removal. These 
estimates show that the yields found in EUV experiments 
with low biases agree with the recomputed yield from the 
SWD experiment. The margin of error is, however, large 
enough that we cannot completely exclude direct influ¬ 
ence from EUV radiation. 


VII. CONCLUSION 

The evolution of an EUV induced hydrogen plasma 
was simulated. The simulations, due to their close cou¬ 
pling to experimental conditions, allowed the magnitude, 
composition, and energy spectrum of the flux from the 
plasma to the surface to be estimated. Our model suc¬ 
cessfully computes the charge bias characteristics of the 
sample to an accuracy of a factor of two. This was used 
to quantitatively compare carbon etch rates between dif¬ 
ferent EUV plasma conditions, as well as compare car¬ 
bon etching under EUV-induced plasma to etching under 
SWD plasmas. In addition, the model describes the fo¬ 
cusing effect of the dielectric surround, which allows the 
measured carbon etch profiles to be better understood in 
a quantitative manner. 

It was observed that the carbon etching mechanism at 
low bias and pressure was different than that at high bias 
and pressure. For the higher energy range, the carbon re- 
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moval yield in EUV-induced plasma was larger than one 
carbon atom per ion. Most probably, etching is domi¬ 
nated by a reactive ion etching process, which may be 
due to the production of methane radicals (e.g. CH^,) 
that can desorb under energetic ion flux. 

However, at low bias energies, the etch rates found 
in EUV-induced plasma agree with those found in SWD 
plasma to within model and measurement uncertainties. 
Therefore, the role of the EUV radiation in a previous 
study [10] was overestimated at low bias energies and 
pressures. 

However, our study does not exclude EUV enhancing 
etch rates at low biases or having a small (< 0.2 C/ion) 
enhancing effect over the range of biases and pressures 
that we studied. Therefore, experiments to determine the 
effect of EUV in the limit of low-to-no bias are required. 
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Figure 8: Experimentally observed and simulated carbon re¬ 
moval as function of pressure, bias and radius. For the con¬ 
ditions of 11.2 Pa, -200V the carbon layer was completely re¬ 
moved, exposing the substrate, at the center of the EUV spot. 
The etch yield for this case was set to the 1.6 C atom/ion as 
in table I. For the conditions of 60 Pa, -200V, the carbon layer 
was completely removed from entire area exposed to plasma, 
hence the case is not shown. The simulated profiles for -50 V 
and -100 V bias were computed via convolution of the simu¬ 
lated ions fluxes with the yield from SWD experiments Fig. 6. 
The profiles were normed to the ion dose per pulse from ta¬ 
ble I. The EUV dose was 10' pulses per sample in all cases, 
which is approximately 170 J. The EUV radiation was con¬ 
centrated in a spot with radius approximately 2 mm. The 
detailed intensity distribution of the EUV spot is not accu¬ 
rately known. 















































